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Abstract 

The decay width and the mass of the D s i(2536) ± have been measured via the decay channel 
Df x -> D*±K% using 232 fb -1 of data collected with the BABAR detector at the PEP-II asymmetric- 
energy e + e~ storage ring. The result for the decay width is r(Z^|) = (1.03 ± 0.05 ± 0.12) MeV/c 2 , 
with the first error denoting the statistical uncertainty and the second one the systematic uncer- 
tainty. For the mass, a value of m(D%) = (2534.85 ± 0.02 ± 0.40) MeV/c 2 has been obtained. The 
systematic error is dominated by the uncertainty on the D*^ mass. The mass difference between 
the Ufi and D*± has been measured to be Am = (524.85 ± 0.02 ± 0.04) MeV/c 2 . 
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1 INTRODUCTION 



The discovery of the mesons D*j(2317) + and -D s j(2460) + [1, 2], with masses considerably lower 
than predicted by potential models [3, 4], has renewed experimental and theoretical interest in the 
spectroscopy of charmed mesons. For a complete understanding of the charmed strange meson 
spectrum, a comprehensive knowledge of the parameters of all known Df mesons is mandatory. In 
this analysis, a precision measurement of the mass and the decay width of the meson D s i(2536) + 
has been performed. The mass is currently reported by the PDG with a precision of 0.6MeV/c 2 , 
while only an upper limit of 2.3MeV/c 2 is given for the decay width [5]. These values are based on 
measurements with 20 times fewer reconstructed candidates compared to this analysis. The 
BABAR experiment, in addition to its excellent tracking and vertexing capabilities, provides a rich 
source of charmed hadrons, enabling an analysis of the with high statistics and small errors. 

Since the uncertainty of the D* + mass is large (0.4MeV/c 2 [5]), we perform a measurement of 
the mass difference defined by 

Am(D+ 1 ) = mCD+i) - m(D* + ) - m{K° s ). (1) 

Additionally, due to the correlation between the masses, the signal in the mass difference 
spectrum is much more narrow than the one from the mass spectrum alone. 

2 THE BABAR DETECTOR AND DATASET 

The data sample used in this analysis corresponds to 232 fb _1 collected with the BABAR detector at 
the PEP-II storage ring from e + e~ collisions at or just below the T(4S) resonance. Furthermore, 
1.16 million DJ[ Monte Carlo events were generated for each of the two decay modes which are 
used for the determination of the detector resolution model. Finally, for resolution studies, D° and 
Kg samples were analyzed using 5% of the main data sample and 20 million simulated e + e~ — > cc 
generic Monte Carlo continuum events. 

The BABAR detector is described elsewhere [6] in detail. Charged particles are detected with a 
combination of a five-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber 
(DCH) filled with a mixture of helium and isobutane, both embedded in the 1.5 T solenoidal 
magnetic field. The transverse momentum resolution is approximately <J Pt /pt = 0.0013(pt/GeV/c)(B 
0.0045. Charged particle identification is done via energy loss measurement within the SVT and 
DCH, and via the Cherenkov light detected in a ring imaging Cherenkov detector (DIRC). Photons 
are detected with a CsI(Tl) electromagnetic calorimeter consisting of 6580 CsI(Tl) crystals. The 
instrumented flux return (IFR) contains resistive plate chambers for the identification of muons and 
long-lived neutral hadrons. For the event simulation we use the Monte Carlo generator EVTGEN [7] 
with a full detector simulation that uses GEANT4 [8]. 

The most critical aspect of this analysis is the quality of the track reconstruction. The track- 
finding algorithm is based on tracks found by the trigger system and by standalone track recon- 
struction in the SVT and DCH. The parameters of a given track are determined using a Kalman 
filter algorithm [9], which makes optimal use of the hit information and corrects for energy loss 
and multiple scattering in the material traversed and for inhomogeneities in the magnetic field. 
The material-traversal corrections change the track momentum according to the expected average 
energy loss and increase the covariance for track parameters to account for both multiple scattering 
and the variance in the energy loss. The latter depends on the particle velocity, so each track fit is 
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performed separately for five particle hypotheses: electron, muon, pion, kaon and proton. A sim- 
plified model of the BABAR detector material distribution is used in the Kalman filter algorithm; a 
wrong simulation of the material will result in a wrong energy loss correction in the reconstruction. 
Besides a good description of the tracking region, a detailed knowledge of the magnetic field is 
also essential for a precise track reconstruction. The solenoid field itself has been well measured. 
The field strength in the tracking volume is known to an accuracy of 0.2 mT. More uncertain is 
the contribution of the solenoid field to the magnetization of the permanent magnets used for final 
focusing and bending of the beams. This is discussed in more detail in section 4.3. 

3 ANALYSIS METHOD 

3.1 D+ candidate reconstruction 

Two decay modes are reconstructed: in both modes the decays to D* + K®, with the K® decaying 
into n + ir~ and the D* + into D°ir + . The D° decays either into K~ir + or K~-k + 7t + it~ . In the following, 
we refer to these two decay modes as K4ir and KGir, respectively, where also the charge conjugated 
states are included. 

For the D° candidate reconstruction, charged kaon and pion candidates are combined to form 
K~ir + (K~ir + Tr + iT~) final states. Fits to the D° mass spectra yield a mean value of 1863.6 MeV/c 2 
and a signal resolution of 7.4MeV/c 2 (8.1MeV/c 2 ) for the K4ir (K6tt) decay mode. The signal 
region is chosen as a mass window of ±18 MeV/ c 2 (±14 MeV/ c 2 ) centered around the obtained 
mean value. The D° candidates are each combined with an additional charged pion to form D* + 
candidates. The fits to the D* + -D° mass difference yield a mean value of 145.4 MeV/c 2 and a signal 
resolution of 0.19MeV/c 2 (0.24MeV/c 2 ) for K4n (K6ir). The signal region is chosen as a mass 
window of ±1.5 MeV/ c 2 centered around the derived mean value for both decay modes. Finally, K® 
candidates are created from oppositely charged tracks. From the fits to the K® mass distributions 
one obtains a mean value of 497.2 MeV/ c 2 and a resolution of 2.5 MeV/ c 2 for both decay modes. For 
the further analysis we define the K® signal region as a ±10 MeV/ c 2 mass window centered around 
the signal mean value for both modes. 

The background of the K® spectrum is further reduced by restricting the angle between the 
Kg direction of flight and the line connecting the primary vertex and the vertex to values 
smaller than 0.15 rad. The D* + candidates are finally combined with the K® candidates to form 
candidates. In order to suppress combinatorial background, a momentum p* > 2.7GeV/c in 
the center of mass system (CMS) is required for candidates. In addition, this restricts the 
source of the candidates to e + e~ — ► cc continuum production. A kinematic fit is applied to 
the candidates which satisfy the above selection criteria in such a way, that the tracks for 
each composed particle have the same origin. The position of the Df x vertex is required to be 
consistent with the e + e~ interaction region. Since the mass difference Am(Df 1 ) is measured, no 
mass constraint is applied. The probability of the vertex fit is required to be greater than 0.1%. 
Initially, more than one candidate per event is reconstructed. Although the multiple use of 
tracks within one reconstructed decay tree is excluded, there might be multiple candidates 
sharing the same daughter candidate. After applying all selection criteria, the average multiplicity 
of candidates per event is 1.008 and 1.02 for the two decay modes. The selection efficiency is 
16% for the K4ir decay mode and 11% for the K6n mode. 

The resulting mass difference spectra Am(Df 1 ) for MC and data are shown in Section 3.3 
and 3.4, respectively. A double Gaussian is fitted to the data spectrum as a rough estimate of the 
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width of the signal. The results for the total width, calculated by adding the weighted Gaussian 
widths in quadrature, are 1.55 ± 0.15MeV/c 2 and 1.66 ± 0.26MeV/c 2 for mode K4n and K6n, 
respectively. Note that for this preliminary fit the intrinsic width and the resolution have not been 
taken separately into account. 

3.2 Resolution model 

Although very clean signals with more than 2400 (2900) entries have been obtained for decay 
mode KAir (K6tt), it is not feasible to obtain the resolution model and the intrinsic width from 
a single fit to data with all parameters allowed to vary. Instead a resolution model is derived 
from the corresponding Monte Carlo samples. The generation of e + e~ — > cc fragmentation 
events with high accuracy is a difficult task, so deviations between simulated and real data are 
possible. In particular the distribution of the CMS momentum p* of the reconstructed Df x mesons 
differs between data and Monte Carlo events. In order to extract a reliable resolution model from 
Monte Carlo, the model is determined as a function of p*. Another method to compensate for 
the inaccuracies of the simulation is to weight the Monte Carlo spectrum according to the p* 
distribution obtained from data and extract the resolution model from the weighted Monte Carlo 
sample. Since the second method relies on both real data and Monte Carlo data and is sensitive 
to the chosen p* binning, it is used as a systematic check and the resolution model is derived from 
the first method. 

The resolution can be extracted by calculating the difference of the invariant mass of a recon- 
structed candidate and the corresponding generated mass. The derived distribution consists only 
of the deconvolved resolution part of the signal. Since for the measurement of the mass the 
mass difference Am(Df 1 ) is taken, the mass difference Am g (D^ 1 ) of the corresponding generated 
candidates has to be subtracted to obtain the deconvolved distribution: 

Am res = Am(D+) - Am g (D+). (2) 

The procedure is the same for both decay modes. The Monte Carlo sample is divided into 25 p* 
bins with a width of 0.07 GeV/ c in the range from 2.70 GeV/ c to 4.45 GeV/c. An unbinned maximum 
likelihood fit of Am res is performed for each bin of p* using a probability density function (PDF) 
assembled from Gaussian functions with fit parameters Am reS] o, r and a®: 

R{Am res ) = / — e 57> d a, (3) 

J <to TO 

where Am reS) o is the mean value and the width is integrated from a minimum value of oq up to a 
maximum of voq. The scale parameter for the upper width limit r is determined for each p* bin, 
but does not vary drastically with p* . Fixed values for r are obtained from a fit with a constant 
function to the r distribution which yields r = 5.64 ± 0.05 (6.40 ± 0.06) for decay mode KAir 
(K6tt). A non-constant 1 st order polynomial can be fitted to the r distribution for mode KAtt; this 
scenario is investigated in Section 4.2. The fits to Am res are repeated with r fixed to the constant 
value obtained which leaves oq as the only free parameter for the width. For each p* bin, oq is 
recalculated. The new oq distribution can be best parameterized by the second order polynomial 
below with coefficients hi. 

MP*) = b + hp* + b 2 p* 2 . (4) 



10 



3.3 Validation of the resolution model 

To verify that the resolution model is reliable, we applied it to the Am(-D^) distribution of the 
Monte Carlo data sample 

Am(D+) = S(A»(B+),T(Dti)) * R(Am res ). (5) 

The generated mass difference distribution S(Afi(D^ 1 ), F(Df 1 )) in the simulation follows a non- 
relativistic Breit-Wigner distribution with mean value An(D+) and width T(L>+) (Table 1). The 
fit of this Breit-Wigner convoluted with the resolution function to Monte Carlo data must return the 
generated values for the mass difference A//(D^) and the decay width T(Df 1 ). Since the convolution 
of the resolution model and the Breit-Wigner function cannot be handled analytically one has to 
use numerical integration methods instead. To compute the convolution integral the Trapezoid Sum 
Rule method has been applied. The convolution window has been chosen as ±10MeV/c 2 , which is 
about 10 times the width of the resolution function, and was divided into 200 bins. The applied 
method returns stable fit results for reconstructed Monte Carlo and reproduces the input values 
for the mass difference and width in the simulation with sufficient accuracy (Fig. 1, Table 1). The 
differences between the reconstructed values and the generated values are assigned as systematic 
uncertainties, which are — 7keV/c 2 (— 14keV/c 2 ) for the mass difference and — 2keV/c 2 (— 9keV/c 2 ) 
for the width . 

Table 1: Results of the test fits to the MC data (statistical errors only) for both decay modes, 
compared with the generated values. 



Parameter 


K4tt 


K6ir 


generated 


Afi{D+) /MeV/c 2 


27.737 ± 0.003 


27.730 ± 0.003 


27.744 


T^/MeV/c 2 


0.998 ± 0.005 


0.991 ± 0.007 


1.000 




Figure 1: Fit of the convolution of the non-relativistic Breit-Wigner function and the resolution 
function to the Am(Df 1 ) distribution in MC (dots) as a crosscheck for the p*-dependent resolution 
model. Left: decay mode KAtt; right: decay mode K6n 
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3.4 Fit to the data 



The assumption that the D~£[ lineshape S(Afi(Df 1 ),T(Df 1 )) follows a non-relativistic Breit-Wigner, 
as used for the MC, is not sufficient for the measurement of the intrinsic width. A better and com- 
monly used description of a resonance lineshape is the following, relativistic Breit-Wigner function: 

[m z — nriQ) + mgl z 

To measure the mass difference A^(D^ L ) and the intrinsic width T(Df 1 ), the extracted p*- 
dependent resolution model R(Am res ) is convoluted with the relativistic Breit-Wigner distribution 
and then fitted to data using an unbinned maximum likelihood fit. For the convolution the same 
numerical integration method as described for the MC test fit (Section 3.3) is applied. The back- 
ground is described by a linear function. The fit parameters obtained for both decay modes are 
listed in table 2. The corresponding mass distributions are shown with the fits in Fig. 2. 



Table 2: Results of the fit to the data (statistical errors only) for both decay modes. 



Parameter 


K4ir 


K6ir 


Afi(D+) /MeV/c 2 


27.209 ± 0.028 


27.180 ±0.023 


r( J D+)/MeV/c 2 


1.112 ±0.068 


0.990 ± 0.059 


Signal yield (events) 


2401 ± 47 


2959 ± 51 




Figure 2: Fit to the data (dots) of the relativistic Breit-Wigner function convoluted with the p*- 
dependent resolution function to the Am(Df 1 ) spectrum to obtain the mass difference A/i(Z)J L ) 
and width T(Df 1 ). The background is described by a first order polynomial, shown by the dotted 
line. Left: decay mode K4ir; right: decay mode KQtt. 
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4 SYSTEMATIC UNCERTAINTIES 



Sources of systematic uncertainties can be divided into three categories: a) discrepancies between 
the resolution in data and Monte Carlo (Section 4.1); b) extraction of the resolution model and 
the fit procedure (Section 4.2); c) uncertainties due to inaccuracies in the track reconstruction 
(Section 4.3). The results for the different systematic errors are listed in Table 3. 

4.1 Resolution studies 

A D° sample including both decays D° — ► K~tt + and D° — > K~ir + 7r + 7t~ and a K® — > it + it~ sample 
have been extracted from data and generic cc Monte Carlo and are used for resolution studies. The 
D° and K® mesons have been reconstructed in the same way as described in Section 3.1, but without 
a requirement for the origin of the particles. Due to the negligible width of these two mesons, the 
resolution can be directly obtained from a fit to the respective signals. To avoid momentum- 
related discrepancies, the samples are divided into p*-dependent subsets. The resolution function 
from Eq. 3 is fitted to each data and cc Monte Carlo subset. The full width at half maximum 
(FWHM) is calculated for each bin as a measure of the width of the resolution function. For the 
ratio between the MC and the data resolution, values of 0.90 ±0.01 for the D° as well as for the 
samples are obtained by fitting a constant function to the p*-dependent ratio. In order to measure 
the impact of this deviation on the results for the Df x signal lineshape in data, the width of the 
resolution function obtained in Section 3.2 is enlarged by 10%. The additional resolution is assumed 
to follow a Gaussian distribution. To achieve this widening, the Am res spectrum (divided into p* 
bins as used for the standard fit explained in Section 3.2) is smeared by a Gaussian distribution 
with a certain width a e which yields a spectrum with a FWHM 10% larger than the original one. 
The parameters of the resolution function are obtained as described in Section 3.2. With this 
resolution function, the fit to data is repeated. As expected, the mass is not affected while the 
width is lowered by 67keV/c 2 (mode K4ir) and 52keV/c 2 (mode KQtt). 

4.2 Fit procedure 

Systematic uncertainties concerning the applied method to extract the mass and width have 
been studied. 

Resolution model As described in Section 3.3, the resolution model has been tested on the MC 
data sample. The small differences between the reconstructed values for Afi(Df 1 ) and T(Df 1 ) and 
the respective generated ones are assigned as systematic errors. 

In a first check, the parameter r of the resolution function is alternatively parameterized by a 
I s * order polynomial. With r now depending on p* rather than being constant, the parameter do 
is estimated as described in Section 3.2. Using the new resolution model, one obtains a negligible 
shift of the mass and — 25keV/c 2 (+7keV/c 2 ) for the decay width. 

Another check involves varying the parameter r within one standard deviation Sr of its statistical 
error. Parameterizing ao using the modified r value, the fit to data is repeated with these resolution 
parameters. As a conservative estimate the larger value of the relative change observed for r — 5r 
and r + 5r is taken as an uncertainty. The mass is not affected while the width changes by — 9keV/c 2 
(K4ir) and -27keV/c 2 (K6it). 
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Numerical integration for convolution Lowering and enlarging the integration range by 
±lMeV/c 2 while keeping the bin width constant does not yield significant changes. Doubling 
and halving the bin size while using the standard 20 MeV/ c 2 integration interval yields no relative 
change of A/x(D^) and T(D S \) either, compared with the standard values. 

Background parameterization In the fit to the Am(D^) distribution for data the combinato- 
rial background was parameterized using a first order polynomial. As a systematic check, the latter 
is replaced by a second order polynomial which yields a mass shift of — 3keV/c 2 (< 0.5keV/c 2 ) and a 
modification to the width by — 7keV/c 2 (— 27keV/c 2 ). In a second test, the background was param- 
eterized by a power-law distribution, e.g. using BG(Am(Df 1 )) = ao + aiAm(Df 1 ) + a,2Am(D~£~ 1 ) a3 . 
Both the mass difference and the width do not noticeably change, compared with the standard 
values. Thus, the deviations from the fit with the 2 nd order polynomial are taken as the systematic 
uncertainty arising from the background parameterization. 

Mass window Enlarging the boundaries of the fitted Am(D^ 1 ) standard region 

(0.015 — 0.045 MeV/ c 2 ) by 10 MeV/ c 2 has no significant impact on the A/x(D^) value and changes 

the width slightly by +4keV/c 2 (-13keV/c 2 ). 

p* correction Two strategies - the parameterization and the weighting method - have been 
followed to take the observed discrepancy between the p* distributions for data and MC data into 
account. The results obtained from both approaches are consistent within their statistical errors, 
giving confidence in the results obtained from the parameterization method. The deviation between 
the two methods can be assigned as the systematic uncertainty arising from the applied correction 
method, whereby only the width is noticeably affected by — 21keV/c 2 (— 9keV/c 2 ). 

4.3 Detector conditions 

In the third category of systematic uncertainties, sources arising from charged particle track recon- 
struction have been studied. The amount of material located inside the inner tracking volume (SVT 
and DCH) as well as the understanding of the B field inside that volume plays an important role 
for the reconstruction of charged particles. In addition to this, a dependency of the reconstructed 
mass difference on the p* momentum and on the azimuthal and polar angles 9 and (f) of the 
has been studied. 

Material inside tracking volume It is possible that the amount of material traversed by 
charged particles in the tracking volume is underestimated [10]. This would result in an incorrect 
energy loss correction for the tracks. Two scenarios have been investigated. First, it is assumed 
that the amount of material inside the total volume (SVT and DCH) is underestimated by 10%. 
Second, it is assumed that only the amount of SVT material is underestimated, by 20%. To study 
the effect, the density of the material at the corresponding detector component is increased by the 
respective number. With the new detector conditions, the data are reanalyzed and the fit to data 
is repeated using the standard resolution model with the parameters obtained in Section 3.2. As 
a conservative estimate, the largest deviations of the mass and width obtained from these fits are 
taken as the systematic uncertainty arising from the tracking region material modification. The 
results are shown in Table 3. 
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Magnetic Field The main component of the magnetic field, which is an essential component for 
the momentum measurement, is the solenoid field which is itself well measured. More uncertain is 
the contribution of the solenoid field to the magnetization of the permanent magnets used for final 
focusing and bending of the beams. The magnitude of the solenoid B field is changed by ±0.02% 
for the systematic study. The magnetization is varied by ±20% in order to account for differences 
between the direct field measurements and the permeability measurements. The largest deviations 
observed for the rescaled solenoid field and the rescaled magnetization are added in quadrature and 
used as a conservative estimate of the systematic uncertainty arising from the B field. 

Length scale Another source of uncertainty for the momentum arises from the distance scale. 
The position of the wires in the DCH is known with a precision of 40 fim. With a drift chamber 
radius of ~ 40 cm, this yields a relative precision of 0.01%. As an estimate for the uncertainty of the 
momentum due to the length scale, a systematic error half the size of the uncertainty obtained from 
the ±0.02% variation of the solenoid field is assigned. For the mass difference, this yields a shift of 
-5keV/c 2 (-8keV/c 2 ) for mode K4it (K6it). The width is shifted by ±lkeV/c 2 (±6keV/c 2 ). 

SVT alignment Detector misalignment can affect the measurement of the angles between the 
tracks and possibly the track momenta. This effect has not yet been studied in this analysis. As 
a conservative estimate, the results obtained from SVT misalignment studies performed for the A c 
mass measurement [10] will be used as a systematic uncertainty for the mass. This yields an 
additional error of ±23keV/c 2 . 

Track quality For the standard Df x candidate selection as described in Section 3.1, no lower 
limit has been set for the number of drift chamber hits left by a charged particle. To estimate the 
effect of improved track quality, the selection of the candidates has been repeated, requiring 
that only those 7r and K candidates are used for the reconstruction which have left at least 
20 hits in the drift chamber. The determination of the resolution model has been repeated using 
the new data samples. The mass difference is shifted by ±5keV/c 2 (— 14keV/c 2 ), while the width 
is shifted by — 105keV/c 2 (— 92keV/c 2 ) with respect to the standard values obtained in Section 3.4. 
In a second test, the tt originating from the D* + decay has been excluded from the tighter track 
selection. Because of its small transverse momentum, it will not traverse large parts of the drift 
chamber. The mass difference is shifted by — 8keV/c 2 (< 0.5keV/c 2 ), while the width is shifted 
by — 110keV/c 2 (— 57keV/c 2 ). As a conservative estimate, the deviations obtained from the first 
measurement will be used as a systematic uncertainty. 

Charge dependence The mass difference A/u,(Df 1 ) and the width T(D^ 1 ) have been determined 
separately for both charges of the D s \. The values obtained are consistent within errors with the 
results for the complete data sample (Table 4). Comparing the different values obtained with the 
respective standard values yields for the mass difference a \ 2 of 1.31 (0.34) for mode K4ir (K6tt) 
and for the width a X 2 of 2.08 (0.85). 

Momentum dependence The mass difference and the width is measured for five different p* 
bins in the range between 2.7 and 4.7GeV/c. The results obtained do not vary drastically compared 
with the standard mass difference value obtained from the complete data sample. The weighted 
mean values for the mass difference and the width lie within the error range of the respective values 
for the full data sample (Table 4). Comparing the p* distributions with the respective standard 
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values yields for the mass difference a x 2 °f 0.76 (3.09) for mode K4ir (K6tt) and for the width a 
X 2 of 3.11 (0.36). 

Angular dependence A dependency of the reconstructed mass on the azimuthal angle <fi has 
been observed using the high statistic D data sample created for the resolution studies, while 
no such effect is visible in the corresponding cc MC. The <\> dependency of the reconstructed D° 
mass follows a sine wave and averages to zero when running over the complete data and thus over 
the whole <fi range. Furthermore, the mass difference A / u(D s l 1 ) and the width r(D s l 1 ) have been 
measured in bins of the angle 9 of the and in bins of the azimuthal angle (f>. In both cases, 
the values obtained for the different bins are consistent with the values obtained from the fits to 
the complete data sample. The weighted mean values all lie within the error range of the standard 
values for Afi(Df 1 ) and T(Df 1 ), respectively (Table 4). Comparing the 9 distributions with the 
respective standard values yields for the mass difference a x 2 °f 5.54 (0.97) for mode A47T (K6tt) 
and for the width a x 2 °f 0-98 (6.76). For the 4> distributions one obtains for the mass difference 
a x 2 °f 5.41 (2.87) and for the width a x 2 of 3.28 (0.37). Larger x 2 values are due to bins with 
greater deviations, caused by the small number of candidates available in this range. 

Run dependence Finally, the mass difference and the width have been determined separately 
for different data taking periods. The results obtained for A/x(D^) and T(Df 1 ) lie within the 
error range of the standard values for the mass difference and the width, respectively (Table 4). 
Comparing the different values obtained with the respective standard values yields for the mass 
difference a x 2 of 0.55 (3.36) for mode K4tt (A67r) and for the width a x 2 of 1.16 (0.39). 

In summary, for both the mass difference A/x(D^) and the decay width T(Df 1 ), no significant 
dependency on the momentum, charge and angle or on the time of data-taking have been 
observed. The fits to the different data subsamples return values within the error range of the 
results obtained from the fits to the complete data sample. 

4.4 Techniques for combining results 

The results obtained for the two decay modes are combined using an expanded x 2 method that 
allows different correlations between the individual components of the systematic uncertainties. As 
a crosscheck, the results are combined using a Best Linear Unbiased Estimate (BLUE) technique 
[11] which takes correlated and uncorrelated errors separately into account. The results obtained 
from the methods are consistent with each other. Since the BLUE method delivers only one single 
combined error, the results from the first method are reported. 

5 SUMMARY 

We have presented a high precision measurement of the mass and the decay width of the meson 
Z) s i(2536) + using the decay mode — > D* + K®. The mass difference between and D* + K® 
for the two reconstructed decay modes is measured to be 

A^D+^k^ = 27.209 ± 0.028 ± 0.031 MeV/c 2 , 
A/j,(D+) K&K = 27.180 ± 0.023 ± 0.043 MeV/c 2 , 
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Table 3: Summary of systematic uncertainties. 





A/x(L£) / keV/c^ 


Y{D^) 1 keV/c^ 




KAir 


K6vr 


K4vr 


K6lT 


Resolution width +10% 


±1 


< 0.5 


±67 


±52 


MC validation 


±7 


±14 


±2 


±9 


Parameterization of r 


< 0.5 


< 0.5 


±25 


±7 


Variation of r 


< 0.5 


< 0.5 


±9 


±27 


Numerical integration (width) 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


Numerical integration (steps) 


< 0.5 


< 0.5 


< 0.5 


< 0.5 


Background parameterization 


±3 


< 0.5 


±7 


±27 


Mass window 


±1 


< 0.5 


±4 


±13 


p* correction 


< 0.5 


±1 


±21 


±9 


Detector material 


±14 


±24 


±20 


±29 


B field 


±10 


±16 


±7 


±13 


Length scale 


±5 


±8 


±1 


±6 


SVT alignment 


±23 


±23 






Track quality 


±5 


±14 


±105 


±92 


Quadratic sum 


±31 


±43 


±131 


±119 



Table 4: Measurements of Aa(DJi) and T(DJi) in dependence of several variables. The first 
row shows the results from the fits to the full data sample, followed by the charge dependent 
measurements. Lines 4 to 6 list the weighted mean values obtained from p*, 9 and <\> depending 
measurements. The last two rows contain the results for different data taking periods. 







/ MeV/c 2 


rCEft) / MeV/c^ 


Dependency 


KAlT 


K6it 


K4ir 


K6ir 


Full data 


27.209 ± 0.028 


27.180 ± 0.023 


1.112 ± 0.068 


0.990 ± 0.059 




27.176 ±0.041 


27.193 ±0.033 


1.201 ±0.101 


1.042 ± 0.086 




27.239 ± 0.037 


27.166 ±0.033 


1.007 ±0.092 


0.933 ± 0.082 


P*(D+) 


27.196 ±0.073 


27.204 ± 0.058 


1.222 ±0.188 


0.932 ±0.155 




27.211 ±0.053 


27.181 ±0.045 


1.107 ±0.129 


0.996 ±0.116 




27.206 ± 0.062 


27.179 ±0.053 


1.108 ±0.151 


0.992 ±0.133 


Runl + 2 


27.236 ± 0.041 


27.232 ± 0.036 


1.022 ± 0.098 


0.947 ± 0.092 


Run?, ± 4 


27.188 ±0.039 


27.145 ±0.031 


1.167 ±0.098 


0.958 ± 0.078 
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with the first error denoting the statistical uncertainty and the second one the systematic uncer- 
tainty. These results correspond to a relative error of 0.15% for the mass difference. This lies within 
the range of precision achievable with the BABAR detector: the J/ip mass has been reconstructed 
with a relative error of 0.05% [6]. 

Combining the results, while taking the systematic errors including the uncertainties of the D* + 
mass (±0.4 MeV/ c 2 ) and of the K® mass (±0.022 MeV/c 2 ) into account, yields a final value for the 
mass of 

m(flji) = 2534.85 ± 0.02 ± 0.40MeV/c 2 , 

while the PDG value for the mass is given as 2535.35±0.34±0.50 MeV/ c 2 . The error on the measured 

mass is dominated by the uncertainty of the D* + mass. The mass difference between the 
and the D* + follows from these results as 

Am = m{D^) - m{D* + ) = 524.85 ± 0.02 ± 0.04 MeV/ c 2 . 

The decay width is measured to be 

r(-D+ 1 )x47r = 1.112 ±0.068 ±0.131 MeV/ c 2 , 
T{D+-y) Ken = 0.990 ± 0.059 ± 0.119 MeV/ c 2 . 

The final combined value for decay width is 

r(L>+) = 1.03 ±0.05 ±0.12 MeV/ c 2 . 

The result for the mass difference Am = m(Df l ) — m(D* + ) represents an improvement in pre- 
cision by a factor of 14 compared with the current PDG value of 525.3 ± 0.6 ± 0.1 MeV/ c 2 . 
Our result deviates by la from the larger PDG value. The precision achieved is comparable 
with other recent high precision analyses performed at BABAR like the A c mass measurement 
(m(yl c ) = 2286.46 ± 0.04 ± 0.14 MeV/ c 2 ) [10]. Furthermore, this analysis presents for the first time 
a direct measurement of the decay width with small errors rather than just an upper limit, 
which is currently stated by the PDG as 2.3 MeV/ c 2 . 
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